Lieb's soliton-like excitations in harmonic trap 
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solitons, vortices, and topological excitations 
PACS 05 . 30 . Jp - Boson systems 

PACS 03 . 75 . Kk - Dynamic properties of condensates; collective and hydrodynamic excitations, 
superfluid flow 

Abstract - We study the solitonic Lieb II branch of excitations in one-dimensional Bose gas in 
homogeneous and trapped geometry. Using Bethe ansatz equations we calculate the "number 
of particles" and the "effective mass" of a soliton. The frequency of oscillations in a harmonic 
trap is calculated. It changes continuously from its "soliton-like" value uJh/V^ in the high density 
mean field regime to cjh in the low density Tonks-Girardeau regime with uj^ the frequency of the 
harmonic trapping. The phase jump of the order parameter is calculated. 



Introduction. — Recent development of experimental 
techniques has opened an exciting possibility to work with 
ultracold Bose gases in one-dimensional (ID) conditions, 
for example in a set of elongated optical traps [I] and in 
magnetic traps created by solid state chips [2]. This de- 
velopment permitted to verify the theoretical predictions 
in highly controllable experiments. The theoretical inves- 
tigations of one-dimensional bosons have been started by 
Marvin Girardeau in Ref. [3], where the case of an infi- 
nite repulsion was considered. This case is often called 
as "Tonks-Girardeau" (TG) limit, although Tonks con- 
sidered ID Boltzmann gas [J. The next important step 
was made by Lieb and Liniger in Refs. [5j[6] where they 
obtained an exact solution of the problem of ground state 
properties and energy spectrum of a one-dimensional Bose 
gas with the delta-function repulsive interaction (Lieb- 
Liniger gas). Probably, the most surprising result of the 
paper [S| is the existence, besides the phonon-like branch 
of elementary excitations (Lieb I branch), which presence 
was natural to assume in analogy with 3D case, also of 
the second branch (Lieb II branch). This branch exists 
in a finite interval of the momenta \p\/p < tt and its en- 
ergy vanishes, when the coupling constants tends to zero. 
The meaning of this branch is clarified in the TG-limit, 
where the Lieb II spectrum coincides with that of an ideal 
Fermi gas, and corresponds to excitations of holes. The 
meaning of the second branch in the opposite (mean-field 



"Bogolyubov" ) limit of weakly interacting bosons was ex- 
plained by Kulish, Manakov and Faddeev [7] (see also [8]). 
They have shown, that the energy-momentum dispersion 
relation of an excitation of the second branch in this 
limit coincides with the relation, obtained by Tsuzuki [9] 
for a soliton, described by the Gross-Pitaevskii equation 
(GPE). Recently Sato et al. have shown that the spatial 
profile of the order parameter, defined as a proper ma- 
trix element, also reproduces the GPE soliton profile [10]. 
Hence, Lieb II branch of excitations in the intermediate 
regime is a result of a quite non-trivial crossover between a 
topological soliton and an excitation of a fermion-like hole. 
Investigation of the properties of these unusual objects is, 
in our opinion, an interesting and important problem. In 
this paper we will investigate dynamics of Lieb II branch 
of excitations in a gas confined to ID harmonic trap. We 
will assume that the size of the cloud is large in compari- 
son with the healing length. Then one can safely use the 
Local Density Approximation (LDA) for the dynamics. 

Local Density Approximation. In the LDA, dy- 
namics of an excitation is defined by its dispersion law in 
a uniform gas. The most convenient description of the 
dynamics is in terms of energy of the excitation £{Vfij ex- 
pressed as a function of its velocity V and the chemical 
potential /x. For a smooth external potential U{x), the 
LDA energy can be obtained by exchanging the chemical 
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potential by its local value fj. ^ ^ — U{x). This means 
that in the course of the motion of an excitation in the 
presence of the external field, the energy e{V,fi — U{x)) 
must remain constant |lllll2j . Differentiating with respect 
to time and taking into account that dx/dt = V wc obtain 
l/V{de/dV)^dV/dt- ide/dti)v{dU/dx) = 0, or 



dV 
dt 



— = -N, 



dU_ 

dx 



(1) 



where the parameters, characterizing the excitations, 

have, correspondingly meaning of the effective mass me/ / 
and the number of atoms Ng in the excitations. For the 
excitations of the second branch in the Bose gas these 
quantities are negative, thus \Ng\ is the number of atoms 
expelled at creation of an excitation. Quantities ^ and 
^ can be easily calculated in the Bogolyubov regime, 
where according to [3] the energy of a soliton is eiV,^) = 
2/l(/Lt— mX^'^)'^/^/(3cm^/^), where 2c is the one-dimensional 
coupling constant (sec Eq. ((8)) below). Correspondingly 



2mN,. 



1/2 



(4) 
(5) 



Thus, the GPE soliton moves in an arbitrary external field 
as a particle of the mass 2m [12]. If the gas is trapped 
in a harmonic trap with frequency w^, the frequency of 
oscillations will be w = ujh/V^ [13] . 

In the opposite TG limit the energy of an second-branch 
excitation can be presented as £(y,/i) = /-t — mV'^/2 and 



meff 



-1, 



m 



(6) 
(7) 



corresponding to the "hole-like" nature of the excitation 
in this limit. In this case u — LJh- In this paper we will 
calculate the characteristic parameters We/ / , -^s and the 
oscillation frequency w for intermediate strengths of the 
interaction. 

It is worth to notice, that an excitation, moving in an 
external field, has a finite lifetime due to decay to several 
excitations. This effect has been investigated in [T3] for the 
GPE solitons. The probability of the decay is supressed 
at small enough uih- In the following we will not consider 
this effect. 

Lieb's equations. — In the Lieb-Liniger model the 
Hamiltonian is written as 



In the original paper [S] authors used the system of units 
with h = l,m = 1/2. Calculation of the second branch 
of the spectrum of elementary excitations is reduced to 
dc solution of a linear integral equation for the function 

J{k,q) 



K 



2TTj{k,q) - 2c 



-K 



J{r, q)dr 
+ {r — k)'' 



= TT — 2 tan 



q — k 



(9) 

The limit of integration K defines the one-dimensional 
density p (and the value of the dimensionless parameter 
7 = c/p) indirectly, as an integral of the solution of an 
equation similar ([5]), but without g-dependent term on 
the r.h.s. |S]. The dependence K^-y) and the inverse one 
^{K) can be calculated following the methods of [S]. Once 
such relations are known, the sound velocity u can be cal- 
culated according to m = —2'y'^d{K/^)/d'y. We use matrix 
methods to solve Eq. ([9]) and similar integral equations. 
To do so we discretized the integral which than is written 
as a (r, q) matrix. The inverse matrix is calculated and is 
multiplied by the discrete representation of the r.h.s. of 
Eq. ®. 

The knowledge of J permits to calculate the dependence 
of the energy e on the momentum p in the parametric form 
(here q is understood as a free parameter): 



K 



K 



e = ^i-q^ + 2 I J{k,q)kdk, p 



-K 



-K 



J{k, q) dk. 
(10) 

The resulting excitation spectrum e is a function of /x 
and p, instead of ^ and V. It is possible to calculate 
V{p,^) = {de/dp)^ from equations ([TO]) . To calculate Ng 
and TJief / in these variables one can use the relations 



Ng 



de\ f dp 



(12) 



i<k 



Results and discussion. — Probably, the best way 
to experimentally verify our predictions is to measure the 
frequencies of oscillations $7 in a trap in different regimes. 
Figure [T] shows the dependence of the frequency O on 
the interaction parameter 7. One can see that the fre- 
quency continuously increases with increasing 7 from its 
GPE value uJh/V2 to the ideal Fermi gas value uJh- The 
most rapid change takes place at 7 ~' 3. There are dif- 
ferent ways of measuring the oscillations frequency. At 
moderately small values of 7, when a soliton still con- 
tains a large number of atoms, one can directly observe 
its motion, like in the experiments [15] . Instead at 7 ~ 1 
the number of atoms in a soliton is small, \Ng\ ~ 1, one 
might exploit the confinement induced resonance (CIR) 
[l6] in order to change the value of 7 in the course of an 
experiment. Initially solitons can be created by a phase 
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Fig. 1: (Color online) Solid line, frequency of oscillations 
in units of the frequency Uh of the harmonic oscillator as a 
function of the interaction parameter 7; dashed line, asymp- 
totic value in TG limit; dash-dotted linea, asymptotic value in 
GP limit; short-dashed line, perturbative solution of integral 
equations, Q,/uJh = 1 — I/7 + 



> 



1 - 



10' 



10^ 




- - TG regime, -A^^ = 1 
GP regime, -A'^ oc y'' 



10" 



10^ 



10' 



10^ 



Fig. 2: (Color online) Solid line, number of particles Ns in 
a stationary soliton (1/ = 0) as a function of the interaction 
parameter 7; dashed line TG limit, Eq. ((6]); dash-dotted line, 
GP limit, Eq. Q. 



imprinting method at small 7. Later the value of 7 might 
be increased by using CIR and soliton oscillations can be 
excited by a parametric modulation of the trap frequency. 
One expects to observe a resonance at the frequency of 
modulation 51/2. The resonance can be detected by heat- 
ing of the gas. 

The frequency of small oscillations is Vt/iUh = 

Nsm/mef f ■ It depends both on the number of particles 
in the soliton Ns and on the effective mass rueff- These 
quantities must be taken in the limit of a small velocity 

— )■ 0. The dependence of the number of particles in the 
soHton at rest Ns{V = 0) on the interaction strength is 
shown in Fig. [2j We find that for small 7 (GP regime) 
\Ng\ ^ 1 and the soliton is a macroscopic object, however 
\Ns\ becomes of the order of 1 at 7 '--^ 1. 

At small 7 the number of atoms is well described by 
the GP result, Eq. (U). However, the situation is different 
for "fast" solitons with a vanishing momentum p — > 0, i.e. 
with the velocity approaching the speed of sound V u. 
According to Eq. (j4|) the number of atoms in a soliton 
tends to zero, |A^s| cx {u — V)^/'^. However, the calcula- 
tions show that Ns tends to a finite value &i V ^ u. It 
is not difficult to prove, that this is only possible, if the 
dispersion law of the soliton has the following expansion 
at p ^ 

e{p)Kup+^, (13) 

with m* ~ \meff{p = 0)|. This relation is quite non- 
trivial, because presence of the contradicts the GPE. 
Indeed, according to GPE e — up (x p^/^. The existence 
of the p^ term was established by Periera et al. [T7] and 
Imambekov, Schmidt, and Glazman (see Eqs. (50)-(51) 
in |18j). Such a term exists both for upper and lower 



branches of elementary excitations. The effective mass is 
the same in the absolute value for two branches. Simple 
calculation permits to present the result of [18] as 

(14) 

In the GPE regime 7^1, velocity of sound u cx p^^^ and 
the second term disappears. Taking to account that the 
presence of the term in dispersion does not violate the 
relation Ns ~ meff/2, one finds for 7 ^ [18] 

\Ns{p = 0)H ^7-'/' = 1.187-^/^ . (15) 

In Fig. [3] we test the obtained result by showing the de- 
pendence of \Ns{p = 0)17^/'* on 7. One can see a good 
agreement for the coefficient in Eq. ([T5]) in GP limit. 

In Galilean-invariant systems with a non-diagonal long 
range order, there exists an important relation between 
the momentum of a solition, the number of atoms Ns and a 
topological characteristic of the soliton, namely the phase 
jump ifs- This relation must be valid both in a superfluid 
Bose system [19] and in a superfluid Fermi gas [20] . For a 
Bose gas the relation can be written as 

hpifs ~ p ^ mNgV . (16) 

The phase jump as a function of the momentum p for 
different values of 7 is presented in Fig. SJ 

Notice that in the TG limit large 7 the jump (ps — tt 
and does not depend on V. We could not find a simple 
interpretation of this result. 

The derivation of Eq. (|T6]) implies an existence of a 
well-defined complex order parameter. The validity of the 
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Fig. 3: Solid line, number of particles in the soliton Ns for 
zero momentum p — (or soliton velocity V — u) scaled with 
7^/4 as a function of 7. Dashed line, constant value of 1.18 as 
expected from Eq. ((15} in GP regime. 

equation in a ID system at not small values of 7, where 
the order parameter is destroyed due to the long-distance 
phase fluctuations, is an interesting open question. In this 
connection it would be quite interesting to compare these 
results with the results of direct numerical simulation of 
the quantum dynamics of solitons at moderate 7, analo- 
gous to presented in |21| and with the exact calculations 
of the profile of the order parameter by the method of [10] . 

Conclusions. — To conclude, by using the exact Lieb- 
Liniger theory we investigate physical characteristics of 
the motion of soliton-like Lieb II branch of excitations. 
The frequency of oscillations, number of atoms in the 
soliton and phase jump are calculated. These quantities 
experience a continuous change in the crossover between 
Gross-Pitaevskii and Tonks- Girardeau limits. The exper- 
imental possibility of the verification of the calculations is 
discussed. 
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